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Abstract. We embark in a program of studying the problem of better ap- 
proximating surfaces by triangulations (triangular meshes) by considering the 
approximating triangulations as finite metric spaces and the target smooth 
surface as their Haussdorff-Gromov limit. This allows us to define in a more 
natural way the relevant elements, constants and invariants s.a. principal di- 
rections and principal values, Gaussian and Mean curvature, etc. By a "natural 
way" we mean an intrinsic, discrete, metric definitions as opposed to approx- 
imating or paraphrasing the differentiable notions. In this way we hope to 
circumvent computational errors and, indeed, conceptual ones, that are often 
inherent to the classical, "numerical" approach. In this first study we consider 
the problem of determining the Gaussian curvature of a polyhedral surface, by 
using the embedding curvature in the sense of Wald (and Menger). We present 
two modalities of employing these definitions for the computation of Gaussian 
curvature. 



1. Introduction 

The paramount importance of triangulations of surfaces and their ubiquity in 
various implementations (s.a. in numerous algorithms applied in robot (and com- 
puter) vision, computer graphics and geometric modelling, with a wide range of 
applications from industrial ones, to biomedical engineering to cartography and 
astrography - to number just a few) has hardly to be underlined here. In con- 
sequence, determining the intrinsic proprieties of the surfaces under study, and 
especially computing their Gaussian curvature is essential. However Gaussian cur- 
vature is a notion that is defined for smooth surfaces only, and usually attacked 
with differential tools, tools that - however ingenious and learned - can hardly 
represent good approximations for curvature of PL-surfaces, since they are usually 
just discretizations of formulas developed in the smooth (i.e. of class at least C^) 
case} 

Moreover, since considering triangulations, one is faced with finite graphs, or, in 
many cases (when given just the vertices of the triangulation) only with finite -thus 
discrete - metric spaces. Therefore, the following natural questions arise: (A) Is 
one fully justified in employing discrete metric spaces when evaluating numerical 
invariants of continuous surfaces? and (B) Can one find discrete, metric equivalents 
of the differentiable notions, notions that are intrinsically more apt to describe the 
properties of the finite spaces under investigations? One is further motivated to ask 
the questions above, since the metric method we propose to employ have already 
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successfully been used in the such diverse fields as Geometric Group Theory, Geo- 
metric Topology and Hyperbolic Manifolds, and Geometric Measure Theory. Their 
relevance to Computer Graphics in particular and Applied Mathematics in general 
is made even more poignant by the study of Clouds of Points (see LWZL^ . .MDj) 
and also in applications in Chemistry (see |Ti)- 

We show that the answer to both this questions is affirmative, and we focus our 
investigations mainly on the study of metric equivalents of the Gauss curvature. 
Their role is not restricted to that of being yet another discrete version of Gaussian 
Curvature, but permits us to attach a meaningful notion of curvature to points 
where the surface fails to be smooth, such as cone points and critical lines. Thus 
we can employ curvature in reconstructin not only smooth surface, but also surfaces 
with "folds", "ridges" and "facets". 

This exposition is organized as follows: in Section 2 we concentrate our efforts on 
the theoretical level and study the Lipschitz and Gromov-Hausdorff distances be- 
tween metric spaces, and show that approximating smooth surfaces by nets and 
triangulations is not only permissible, but is, in a way, the natural thing to do, in 
particular we show that every compact surface is the Gromov-Hausdorff limit of 
a sequence of finite graphs.^ In Section 3 we introduce the best candidate for a 
metric (discrete) version of the classical Gauss curvature of smooth surfaces, that 
is the Embedding, or Wald curvature. We study its proprieties and investigate the 
relationship between the Wald and the Gauss curvatures, and show that for smooth 
surfaces they coincide, so that the Wald curvature represents a legitimate discrete 
candidate for approximating the Gaussian curvature of triangulated surfaces. Sec- 
tion 4 is dedicated to developing formulas that allow the computation of Wald 
curvature: first the precise ones, based upon the Cayley-Menger determinants, and 
then we develop (after Robinson) elementary formulas that approximate well the 
Embedding curvature. We conclude with three Appendices. In the first Appendix 
we present three metric analogues for the curvature of curves, namely the Menger, 
Alt and Haantjes curvatures and study their mutual relationship. Furthermore we 
show how to relate to these notions as metric analogues of sectional curvature and 
how to employ them in the evaluation of Gauss curvature of triangulated surfaces. 
Next we present yet another metric analogue of surfaces curvature, based, in this 
case, upon a the modern triangle comparison method, namely the Rinow curva- 
ture. We investigate its proprieties and show (following Kirk (j£,)) that in the 
case under investigation the Rinow and Wald curvatures coincide (and therefore 
Rinow curvature also identifies to the Gauss curvature) . The third and last Appen- 
dix is dedicated to the development of determinant formula for the radius of the 
circumscribed sphere around a tetrahedron, with a view towards applications. 

2. The Haussdorff-Gromov limits 

2.1. Lipschitz Distance. This definition is based upon a very simple^ idea: it 
measures the relative difference between metrics, more precisely it evaluates their 
ratio; i.e.: 

The metric spaces {X, dx), (Y, dy) are close iS 3 f : X ^ Y s.t. '^^J'^fl^'^f' ~ 1, 



^ For the relevance of these notions in the study of classical curvatures convergence, see ]CMS[ . 

m 
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Technically, we give the following: 

Definition 2.1. The map / : {X, dx) — * {Y, dy) is bi-Lipschitz iff 3 c, C > s.t.: 

c-dx{x,y) < drifxjy) <C ■ dx{x,y) . 

Remark 2.2. The same definition applies for two different metrics d\,d2 on the 

same space X . 

Definition 2.3. Given a Lipschitz map f : X ^ Y , we define the dilatation of f 
by: 

dilf= sup — — — . 

x^yex ax{x,y) 

Remark 2.4. The dilatation represents the minimal Lipschitz constant of maps 
between X and Y. 

Remark 2.5. If / is not Lipschitz, then dil f = oo. 

Remark 2.6. (1) / Lipschitz / continuous. 
(2) / bi-Lipschitz => / homeo. on its image. 

Remark 2.7. We have the following results: 

Proposition 2.8. Let f,g:X^Y be Lipschitz maps. Then: 
(C') 9 ° f *5 Lipschitz 
and 

(b) dil {g o f) < dil f ■ dilg 

Proposition 2.9. The set {/ : {X, d) — > M | /Lipschitz} is a vector space. 

Now we can return to our main interest and define the following notion: 

Definition 2.10. Let {X,dx), (F, dy) be metric spaces. Then the Lipschitz dis- 
tance between (X, dx) and (Y, dy) is defined as: 

dLiX,Y) = inf logmax (dil f, dil f~^) 

f hi — Lip. 

Remark 2.11. If ^ f bi-Lipschitz between X and Y, then - remembering Remark 

2.2 we put dL{X,Y) = oo (i.e. di, is not suited for pairs of spaces that are not 

bi-Lipschitz equivalent.) 

having defined the distance between two metric spaces we now can define the 
convergence in this metric in the following natural way: 

Definition 2.12. The sequence of metric spaces convergence to the 

metric space {{X,d)} iff 

\imdL{Xn,X) = 

n 

(In this case we write: {Xm dn)'^0). 

Example 2.13. Let St be a family of regular surfaces, Si = ftiU); where U is an 
open set, U = intU C M; such that the family {ft} of parametrizations is smooth 
(i.e. F :U xR^m^ eC^; where F{{u,v),t) =* ft{u,v)). Then d„{St, 80)^-^1^0. 

^Here and in the sequel "fx" etc. ... stands as a short-hand version of "/(a;)". 
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Remark 2.14. If F is not smooth (only continuous) then we do not necessarely have 
that St jZ^^qS'o . 

We have the following significant theorem: 

Theorem 2.15. The cLl satisfies the following conditions: 

(a) dL > 0; 

(b) dL is symmetric; 

(c) d^ satisfies the triangle inequality; 
Moreover, if X, Y are compact, then: 

(d) dL{X,Y) = X <^Y (i.e. X is isometric to Y); 
that is 

is a metric on the space of isometry classes of compact metric 
spaces 

Remark 2.16. Let us recall the following 

Definition 2.17. {Xn,dn)—^{X,d) = du^ d as a real function; i.e. 

sup \dn{x,y) - d{x,y)\—*Q 
x,yex " 

(where "u" denotes uniform convergence.) 

Then Xn^X ^ Xn~£*X but Xn^X ^ Xn^X. However, for finite spaces 
indeed Xn^X ^ X„-jr^X. 

2.2. Gromov-Hausdorff distance. This is also a distance between compact met- 
ric spaces ((distinguished) up to isometry!). However it gives a weafcer topology (In 
particular: it is always finite (even for pairs of non-homeomorphic spaces.) )^ 
We start by first introducing the classical 

2.2.1. Hausdorff distance. 

Definition 2.18. Let A,BC (X, d). We define the Hausdorff distance between A 
and B as: 

dH{A, B) = inf{r > | A C Ur{B), B C Ur{A)} 

(see Fig. 1); where Ur{A) is the r- neighborhood of A, Ur{A) = UaeA-^''(^)' (here, 
as usual: Br{a) = {x & X \ d{a,x) < r}.) 

Another (equivalent) way of defining the Hausdorff distance is as follows: 

dH{A, B) = maxjsup d{a, B), sup d{b, A)} . 

aeA beB 

(see Fig. 2) 

We have the following 

Proposition 2.19. Let {X,d) be a metric space. Then: 

(a) dH is a semi-metric (on 2^ ). (i.e. A = B ^ dniA, B) = O.j 

(b) dH{A,A)=QyA^X. 

(c) {A = AandB = B) ^ {dH{A, B) = ^ A ^ B) . 

i.e. dn is a metric on the set of closed subsets of X. 

'^The relationship between the Lipschitz and the Hausdorff distances is akin to that between 
the C° and norms in Functional Spaces. 
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X 



Figure 1. 

Notation We put: MiX) = {{A <Z X \ A = A}, (Ih) = 2^ /dn ■ 

Remark 2.20. (1) \i X is compact and if {^„}„>i C X is a sequence of compact 
subsets of X, then: 

(a) An+l ^An^Ar,-^ n„>l ■ 

(b) A„ C A„+i ^ Ar,^ U„>1 An ■ 

(2) For general subsets A„^A e A4{X) , and 

(a) A = {lim„ a„ | a„ e A„ ; n > 1} . 

(b) A = ^n„>i(u::=„^™)- 

(3) If An^A, and if the sets A„ are all convex, then A is convex sets. 

We have the following two important results, which we present without their 
respective (lengthy) proofs: 

Proposition 2.21. X complete A4{X) complete. 
Theorem 2.22. (Blaschke) X compacts Ai{X) compact. 




Figure 2. 
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2.3. The Gromov-HausdorfT Distance. We are now able to define the Gromov- 
Hausdorff distance using the foUowing basic guide-hnes: we want to get the maxi- 
mum distance that satisfies the following two conditions: 

(a) dGH{A,B) < dH{A, B)yA, B c X (i.e. set that are close as subsets of X will 
still be close as abstract metric spaces); 

and 

(b) X isometric to Y ^=> dGH{X,Y) = 0. 

Definition 2.23. Let X, Y be metric spaces. Then the Gromov-Hausdorff distance 
between X and Y is defined by: 

dGH{X,Y) = mfdf,{f{X),g{Y)) 

where the infimum is taken over all the isometric embeddings f : X '-^ Z, g : Y ^ Z 
into some metric space Z. (See Fig. 3). 




Figure 3. 

Remark 2.24. li X = S^, with the spherical metric, and Z = M^, with the Euclidian 
metric, then f{X) ^ X (!) 

Example 2.25. Let Y be an e-net'' in X. Then dGH{X,Y) < e. 
Proof Take Z ^ X ^ X' , Y = Y' . 

□ 

Remark 2.26. It is sufficient to consider embeddings / into the disjoint union of 
the spaces X and Y, XY[Y. 

Remark 2.27. (1) X,y bounded => dGH{X,Y) < oo. 

(2) If diamX,diamY < oo, then dGH{X,Y) > ^\diamX — diamY\'^. 

^ Definition Let (X, d) be a metric space, and let A G X. A is called an e-net iff 
d{x,A) <£,'ix & X. 

7 ^ 
diamX = sup^ y^x d{x, y) 
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However , the straightforward definition of doH may be difficult to implement. 

Therefore wc would like to estimate (compute) cIgh by comparing distances in X 
vs. distances in Y (as done in the cases of uniform and Lipschitz metrics). We 
start by defining a correspondence between metric spaces: X < — > Y, given by 
correspondences a; y between points x G X,y G Y. 

Remark 2.28. A correspondence is not necessarely a function, that is to a single x 

may correspond to several y-'s. 

We shall prove that 

(*) dGH{X,Y) < r 3 a correspondenceX < > Y 

s.t. {x ^ y,x' ^ y') =^ \dx{x,x') - dY{y,y')\ < 2r 

Formally, we have: 

Definition 2.29. Let X,Y denote sets. A correspondence X < — > F is a subset of 
the Cartesian product of X and Y: TZ C X x Y s.t. 

(i) Vx € X,3y€ Y, s.t.{x, y) € Tl; 
and 

(ii) Vy eY,3xe X, s.t.{x, y) e n. 

Example 2.30. Any surjective function f : X ^ Y represents correspondence 

n^{{x,,f{x))}. 

Remark 2.31. 7?. is a correspondence <J=^ 3 Z and 3 f : Z ^ X ,3g : Z ^ Y; f,g 
surjective, s.t. TZ = {{f{z),g{z)) \ z e Z}. 

Definition 2.32. Let 7?. be a correspondence between X and Y, where X, Y are 
metric spaces. We Define the distortion of TZ by: 

disTZ = sup{ \dx{x,x') - rfy(?/,j/')l | {x,y) ,{x' ,y') e TZ} . 

(See(*).) 

Remark 2.33. (1) If 7^ = {{x, f{x))} is a correspondence induced by a surjec- 
tive function f : X ^ Y , then disTZ = dis /, where: 

disf= sup \dY{fa,fb) - dx{a,b)\^ . 
a,bex 

(2) If 7^ = {{f{z), g{z))}, where f : X ^ Z, g :Y ^ Z are surjective functions, 
then: 

disTZ= sup {\dx{fz,fz')~dY{gz,gz')\}. 

z.z'ez 

(3) 7?. = iff 7?. is associated to an isometry. 
We bring, without proof, the following theorem: 

Theorem 2.34. Let X,Y be metric spaces. Then: 

dGH{X,Y) = ^mi{disTZ); 

where the infimum is taken over all the correspondences X < — > Y. 

^ Remember that any correspondence can be expressed in this functional manner. 
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Before bringing the next result (which is very important in determining the 
topology ) we first introduce one more notion: 

Definition 2.35. / : X — > F is called an s-isometry (e > 0), iff 

(i) dis f < e, 
and 

(ii) f{x) is an £-net in Y. 

Remark 2.36. / e-isometry =^ / continuous. 

Corollary 2.37. Let X,Y be metric spaces and let £ > 0. Then: 

(i) dGH{X, F) < e 3 2e-isometry f -.X ^Y. 

(ii) 3 s-isometry f : X ^Y => dGH{X,Y) < 2s. 

Proof (i) Let X ^Y s.t. disU < 2s. 

For any x G X and f{x) G Y, choose y = f{x) s.t. {x, f{x)) G TZ. Then x i— > f{x) 
defines a map f : X ^ Y. Moreover: dil f < dilTZ < s. 
We shall prove that f{X) is a 2e-net in Y. 

Indeed, let x G X and y e F s.t. (a;, y) G TZ. Then d{y, fx) < d{x, x) + dis TZ < 2r, 
thence d{y,f{X)) < 2r. 

□ 

Let / be an 2e-isometry. Define TZ C X xY, TZ= {{x, y) \ d{y, fx) < s}. 
Then, since f{X) is an e-nct it follows that 7?. is a correspondence. 
Then V (x, y), {x' , y') e 7?. wc have: 

I dviy, y')-dx{x, a;') I < I d{fx, fx')-d{x, x') \+d{y, fx)+d{y', fx') < dis f+s+s < 3s . 
=^disTZ<3s dGH{x,y) <3r/2 <2r. 

□ 

The next result is of great importance (in particular so in our context): 

Theorem 2.38. dcH is a (finite) metric on the set of isometry classes of compact 

metric spaces. 

Proof It suffices to prove that dcHiX, Y) = 0=^ X ^Y.^ 

Indeed, let X, Y be compact spaces s.t. doH = 0. Then it follows from the previous 

Corollary (for s = 1/n) that 3 (/„)„>i, fn-X^Y s.t. dis fn-^0. 
let S C X, S , 15*1 = Ho . Using a Cantor-diagonal argument one easily shows that 
3 {fnk)k>\ C (/„)„>! s.t. {fnk)>i Converges in Y,\/x G S. Without restricting the 
generality wc may assume that this happens for {fn)n>i itself. Thus we can define 
a function f : X ^ Y hy putting: f{x) = lim„ fn{x). 

But I d{fnX, fny) - d{x, y) I < dis fn^O d{fx, fy) = lim d{fnX, In other 
words f\g is an isometry. But S = S, therefore this isometry can be extended to 
an isometry / from X to Y. In a analogous manner one shows the existence of an 
isometry f : X ^ Y. 

□ 

Remark 2.39. X„-^X X„^X ^n^X. 

In fact, the following relationship exists between and 
Theorem 2.40. Xn—'X s-nets in Xn—*s-nets in X. 

GH L 
^We shall write: X = Y ii X is isometric to Y . 
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One can formulate this assertion in a more formal manner and it directly (see 
IG+K pg. 73). However we shall proceed in more "delicate" manner, starting with: 

Definition 2.41. Let X,Y be compact metric spaces, and let e,6 > 0. X,Y are 
called e-S- approximations (of each-other) iff: 3 {xi}^i C X, 3 {yi}fLi C Y s.t. 

(1) {xi}fLi is an e-net in X and {yi}fLi is an e-net in Y; 
(ii) \dx(xi,Xj) - d(yi,yj) | < (5 Vi, j G {1, ...,7V}. 

An (e, £)-approximation is called, for short: an e- approximation. 

The relationship between this last definition and the Gromov-Hausdorff distance 
is first revealed in 

Proposition 2.42. Let X, Y be compact metric spaces. Then: 

(1) If Y is a {e, 5) -approximation of X , then dGHiX,Y) <2e + S. 

(2) dGH{X,Y) < e =4> Y is a be -approximation of X . 

Proof (1) Condition (ii) of Def. 2.41. is equivalent to disTZxoYo < where 
^0 = {x^jfLi Ay^}fLl■ But disTZxoYo <S =^ dGHiXo,Yo) < S/2. Now, since 
Xq and Yq are e-nets in X, resp. Y, it follows that dcHiXjXo) < e, dG//(y, Iq) < 
e. From here and from the dGH{Xo,Yo) < S/2 follows, by means of he triangle 
inequality, that c?g//(X, Y) < 2e + S. 

□ 

(2) By Cor. 2.37., there exists a 2e-isometrie f : X ^ Y. Let Xq ~ {xi}fLi be an 
e-net, and let yi = f{xi). 

Then | d{xi,Xj) — d^yi, yj) | < 2e < 5e. Therefore suffice to prove that Yq = {yi}^i 
is a 5e-net in Y. 

Indeed, if y G Y, then, since f(X) is an 2e-net in y, 3x E X s.t. d{y, f{x)). Now, 
since Xq is an e-net in X, 3xi G Xq, s.t. d{x,Xi) < e. 
Therefore: d(y, y,) = d(y, /(x,)) < /(x)) + d(/(x), /(x,)) 

< 2e -I- d{x, Xi) + dis f 

<2e-\-e-\-2e <5e. 

□ 

Remark 2.43. Prop. 2.42. (X„ — 'X) <^ (Ve > 0, Xn is an e-approximation, 

GH 

Vn large enough.) 

More precisely we have the following Proposition: 

Proposition 2.44. Let X, {Xn}j° compact metric spaces. Then: 

Xn—^X <^=^ Ve > 0, 3 a finite e-net S <Z X and 3 a finite e-net Sn C Xn, s.t. 

GH 

Sn'Z'XS and, moreover, \Sn\ = \S\, for large enough n. 

GH 

2 42 

Proof ('^=) If 5*, 5*^ exist, then =^ is an ^-approximation of X — 'X\In. 

GH 

□ 

(=>) Let S be an finite e/2-net in X . 
We construct in X„ corresponding nets Sn (to be more precise, we define: Sn = 
fn{X), where /„ is an e„-approximation, fn : X ^ Xn, en 0.) Then Sn^-tS 

GH 

and, in addition, Sn is an e-net in S (for n large enough). 

□ 

We make the following extremely important Remark: 
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Remark 2.45. Let M{n, k, D) be an n-dimensional manifold, of (sectional, Ricci) 
curvature < fc, and s.t. diamM > D. Then (Ai^dcH) is compact. However, it 
should be noted that this result doesn't hold for curvature < fc. (only for Vol{A4) < 
Vq) and injectivity radius > tq. 

Note With the notations of the precedent Proposition, the distances in 5„ con- 
verge to the distances in S, as Xn^X, therefore The Geometric Proprieties 
of Sn will converge to those of S. Thus we can use the Gromov-Hausdorff 
each and every time The Geometric Proprieties of X„ can be expressed 
in term of a finite number of points, and, by passing to the Umit, auto- 
matically obtain proprieties of X. 

A typical example is that of the intrinsic metric i.e. the metric induced by a length 
structure (i.e. path length) by a metric on a subset (of a given metric space). (See 
Fig. 4 for the classical example of surfaces in R'^.) 




Figure 4. 

On a more formal note, we have the following characterization of intrinsic metrics: 

Theorem 2.46. Let {X,d) be a complete metric space. 

(1) If\/ x,y € X, 3^xy, then d is strictly intrinsic. 

(2) If\/x,y G X and Ve > 0, 3 the e-middle of xy, then d is intrinsic. 

Where we used the following definitions and notations: 

Definition 2.47. (1) Given x,y points in {X,d), the middle (or midpoint) of 
the segment xy (more correctly: 'a midpoint between "x" and "y" ') is 
defined as: ^ 

-xy = z, d{x,z) = d(z,y). 

(2) d is called strictly intrinsic iff the length structure is associated with is 
complete. 

(3) Let d be an intrinsic metric, z is an e-middle (or an e-midpoit) for xy iff: 
I 2d{x, z) — d{x, y) \ < e and | 2d{y, z) — d{x, y)\ < e. 

Remark 2.48. The converse of Thm. 2.46. holds in any metric space, more precisely 
we have: 

Proposition 2.49. If d is an intrinsic metric, then ^xy exists, Vx^y. 
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The following Theorem shows that length spaces are closed in the GH-topology: 

Theorem 2.50. Let {X„} he length spaces and let X be a complete metric space 
s.t. X^^X. 
Then X is a length space. 

Proof We have already presented the idea of the proof: it is sufficient to show that 
for every x,y there exist an e-midpoit (Ve > 0). 

Indeed, let n be such that dcH < jq- Then, from the a preceding result, it follows 
that there exist a correspondence Xn^^^X s.t. disTZ < |. 

Let x,y £ Xn, xSx, ySy. Since X„ is a length space, 3z G X„ s.t. z = |- 

7^ — 

midpoint of XnVn- Consider z € X, z^z. Then: 

111 1 I e 2e 

W^A - -\xy\\ < \\xz\ - -\xy\\ + 2dis'R <- + —<£. 
' 2 ' 2 ' 5 5 

(Here we write \xy\ instead of d{x,y), etc.) 

In a similar manner we show that: \\yz\ — ^\xy\ \ < s; i.e. e-midpoit of xy. 
The next Theorem and its Corollary are of paramount importance: 



□ 



Theorem 2.51. Any compact length space is the GH-limit of a sequence of finite 
graphs. 

Proof Let e, (5 ((5 <C e) small enough, and let 5 be a (5-net in X. 
Let G — {V,E) be the graph with V = S and E = {{x,y)\d{x,y) < e}. we 
shall prove that G is an e-approximation of X, for 6 small enough (i.e. for S < 
^diam{X)). (See Fig. 5.) 




Figure 5. 

But, since S is an £-net both in X and in G, and since doix^y) > dx{x,y), it is 
sufficient to prove that: 

dG{x,y) < dx{x,y) + e. 
Let 7 be the shortest path between x and y, and let xi,...,Xn € 7 s.t. n < 
length{-f)/e (and |a;i,a;i+i)| < e/2). Since WxiByi £ S s.t. \xi,yi\ < S, it follows 
that \yiyi+i\ < \xiXt+i \ + 26 < e. (See Fig. 6) 
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Therefore, (for 5 < e/4) 3 an edge e S G, e = yiVi+i- From this we get the foUowing 
upper bound for dcix, y): 

dG{x,y) < Solj/ij/i+il < SolxiXi+il + 2Sn 

But n < 2length{'^) / e < 2diam{X)/e; therefore: 

dcix, y) < \xy\ + S ^ < \xy\ + e . 

e 

(because S < / Adiam{X)) . 




Figure 6. 

So, for any e > 0, 3G = an e-approximation of X. Then, Gn~^X. 

□ 

Corollary 2.52. Let X he a compact length space. Then X is the Gromov- 
Hausdorff limit of a sequence {G„}„>i of finite graphs, isometrically embedded 
in X . 

Remark 2.53. (1) If Gn^X, Gn = If 3No G N s.t. 

(*) \En\<No, VneN, 

then X is a finite graph. 
(2) If condition (*) is replaced by: 

(**) IK I <7Vo, VtiGN, 

then X will still be always a graph, but not necessarily finite(!) 
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3. The Embedding Curvature 



3.1. Theoretical Setting. This is basically a comparison- curvature (as is the more 
"modern" CAT^^ approach). This is done with quadruples instead of triangles 
(like in the Alexandrov-Topogonov method). It is in a sense a more natural idea, 
since quadruples are classically^^ the "minimal" geometric figures that allow the 
difi'erentiation between metric spaces. This allows for a much more easier and 
rapid development of the theory than the triangle-based comparison. Moreover we 
shall show that the two Theories coincide on those metric space on which both 
can be applied, i.e. metric spaces that are (a) "planar" and (b) "rich enough" i.e. 
contain quadrangles, s.a. classical (PL-smooth) surfaces in M.'^}'^ 

Definition 3.1. Let {M,d) be a metric space, and let Q — {pi,...,p4} C M, 
together with the mutual distances: dij = dji = d{pi,pj); 1 < i, j < 4. The set Q 
together with the set of distances {dij}i<ij<4 is called a metric quadruple. 

Remark 3.2. One can define metric quadruples in slightly more abstract manner, 
without the aid of the ambient space: a metric quadruple being a 4 point metric 
space; i.e. Q = [{pi, ..-,^4}, {dij}), where the distances dij veriiy the axioms for a 



Before we proceed to the next definition, let us introduce the following 

Notation 5^ denotes the complete, simply connected surface of constant curvature 
K, i.e. = if K = 0: 5«- = §^ , if k > 0; and = , if k < 0. Here 

Sk = denotes the Sphere of radius R = and = li^^y—j- stands for the 

Hyperbolic Plane of curvature V— '^i as represented by the Poincare Model of the 
plane disk of radius R= 1/ 

Definition 3.3. The embedding curvature k{Q) of the metric quadruple Q is de- 
fined be the curvature k of S^, into which Q can be isometrically embedded. (See 
Figures 7 and 8 for embeddings of the metric quadruple in 5*^ and H,^, respectively.) 

We can now define the embedding curvature at a point in a natural way by pass- 
ing to the limit (but without neglecting the existence conditions), more precisely: 

Definition 3.4. Let (M, d) be a metric space, and let p G M be an accumulation 
point. Then p is said to have Wald curvature k,w{p) iff 
(i) ^ TV gA/'(p), TV linearis ; 



vature of quadruples, then the very existence of k{Q) is not assured, as it 
is shown by the following 

Counterexample 3.6. The metric quadruple of lengths 



i.e. Cartan-Alcxandrov-Topogonov 
^"'"as illustrated by the time-honored prineiples of Projective Geometry... 

In this sense CAT spaces are more "potent" : they can be employed in studying mathematical 
objects that not (neccessarilly) contain quadrangles, e.g. trees, Cayley graphs, etc.. 
The neighborhood iV of p is called linear iff N is contained in a geodesic. 



metric. 




dl2 = di3 = di4 = 1; (^23 = £^24 = £^34 = 2 
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admits no embedding curvature. 

(2) Even if a quadruple has an embedding curvature, it still may be not unique 
(even if Q is not liniar), indeed, one can study the following examples: 

Example 3.7. (a) The quadruple Q of distances dij = 7r/2, 1 < « < j < 4 
is isometrically embeddable both in Sq = and in Si = S^. 




Figure 8. 
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(b) The quadruple Q of distances c?i3 = (ii4 = d2z — ^24 = tt, di2 = 
^34 — 37r/2 admits exactly two embedding curvatures: ki € (1-5,2) 
and K2 = 3. (See (Ml-) 

However, for "good" metric spaces^** the embedding curvature exists and it is 
unique. And, what is even more relevant for us, this embedding curvature coincides 
with the classical Gaussian curvature. The proof of this result is rather long and 
tedious, therefore we shall present here only a brief sketch of it. (This will prove to 
be somewhat redundant anyhow, in view of the more general results presented in 
the previous section, a fact but we shall emphasize later in our presentation).) 
The Main ingredient for this proof, and for the analysis of yet another another 
approach to curvature (the CAT one) is provided by the following string of propo- 
sitions (which are just generalizations of the well known high-school triangle in- 
equalities): 

Proposition 3.8. Let A(pi,gi,ri) C S^.^ and A(p2,92,7'2) C 5^2, s.t. 
Piqi = ^292, Piri = P2r2 and Z(gi,pi,ri) = Z(g2,P2, ?'2)- 
Then: ki < K2 =^ qiVi > g2?'2- 

Proposition 3.9. Let pi,qi,ri € 5^1, P2i<l2,i"2 G two isometric triples of 
points, s.t. the triple pi,qi,riis not linear. Then: 

-^(91,^1, fi) < A'l2,P2,r2),Z{pi,qi,ri) < Z(p2, 92, 7-2), /(gi, r-i,pi) < ^(92, ^2,^2)- 

Proposition 3.10. LetQi {pi, 51, ri, si}, Q2 — {P2, <l2,T2, S2} be non-linear and 
non-degenerate quadruples inS^^, 5^2, respectively. // A(pi, gi, ri) = A(p2,92,?'2) 
and Ki < K2, then: 

(1) pisi = P2S2, qisi = 9252 risi > r2S2 ; 

(2) riSi = r2S2, giSi = g2S2 =^ PiSi > P2S2 ; 

(3) piSi = P2S2, riSi = r2S2 =^ qiSi < 92^2 ■ 

In order that we fully exploit the results above we need the following definition: 

Definition 3.11. A metric quadruple Q = Q(pi,P2,P3,P4), of distances 

dij = dist{pi,pj), i = 1, ...,4, is called semi-dependent (a sd-quad, for brevity), iff 3 

of its points are on a common geodesic, i.e. there exist 3 indices - e.g. 1,2,3 - s.t.: 

^12 + 0^23 = 0^13 • 

Now we can easily formulate the following immediate consequence of Prop. 3.10. : 

Corollary 3.12. A sd-quad admits at most one embedding curvature. 

Unfortunately - as we have already noticed - in the general case the uniqueness 
of the embedding curvature is not guaranteed. However we can be a bit more 
explicit using the following definition: 

Definition 3.13. Let Q = {p, q, r, s} be a non-linear and non-degenerate quadru- 
ple. Q is called planar iff Z{q,p, r) -\- Z{q,p, s) + Z(s,p, r) ~ 2tt. 

Then we have 

Proposition 3.14. Let Q — {p,q,r,s} be a a non-linear and non-degenerate 
quadruple in iS^. Then 

(1) Lf Q is planar, then it admits no isometric embedding in ni > k- 



i.e. spaces that are locally "plane like' 
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(2) If Q is not planar, then it admits no isometric embedding in 3^,2, ^2 < k- 

Corollary 3.15. Let Q = {p, q, r, s} be a a non-linear and non- degenerate quadru- 
ple. Then Q has at most two different embedding curvatures. 

In fact we can state a much stronger assertion, of which Example 3.7. (a) is just 
a very particular case: 

Proposition 3.16. Vp S S^, and Vk > 0, 3U G Af{p) s.t. 3 a nonlinear, 
non-degenerate quadruple Q d U of embedding curvature 0. 

Proof. Let 71,72 G U, two great-circle arcs s.t. 71 H 72 = p. 
Let gi, 92 G 71 s.t. pqi — pq2 ^ and let g G 72 s.t. pq < 7r/2^/K. 
Consider l^{q[q'2,q') C R^, l^{q[q^2,q') = A{qiq2,q), let p' = y.q!^, and let 
h = q'p'. 





Figure 9. 

Then since < k, Proposition 3.10.(3) applied to the quadruples {q,qi,q2,p} 
and {q' , q'l, q'2,p'} implies that h < pq. 

Now let X £ -f2, X between p and q, and let x' € s.t. A{q[q2,x') = A{qiq2,x) 
s.t. X and q' are on different sides of the line qiq2- Then, x ~ p ^ xq > x'q' , 
and X = q ^ xq = < x'q' = 2h, where, in this case x' ~ q". (See Figure 9.) 



i.e a quarter of the length of a great circle in 
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Then it follows from a continuity argument that 3 G 72, a^o between p and s.t. 
xoq = x'^q', thus implying that {gi, 92, g, a;} = {gi , 92^ x'}. □ 

Remark 3.17. {51,(72,9,2:} is planar, while {gj^, 2:'} is not planar. 

3.2. The Wald Curvature vs. Gauss Curvature. The discussion above would 
be nothing more than a nice intellectual exercise where it not for the fact that the 
metric (Wald) and the classical (Gauss) curvatures coincide whenever the second 
notion makes sense, that is for smooth (i.e. of class > C^) surfaces in . More 
precisely the following theorem holds: 

Theorem 3.18. (Wald) Let 5 C K^, S GC"',m>2 be a smooth surface. 
Then Kwip) exists, for all p G S, and nwip) = '«g(p),Vp S M. 

Moreover, Wald also proved that a partial reciprocal theorem holds, more pre- 
cisely he proved the following: 

Theorem 3.19. Let M be a compact and convex metric space. 

LfKwip) exists, for allpG M, then M is a smooth surface and kw{p) = i^g{p),^P^ M. 

Remark 3.20. I one tries to restrict oneself, in the building of Definition 3.4. only 
to sd-quads, then Theorem 3.19. holds only if the following presumption is added: 

Condition 3.21. AI is locally homeomorphic to M.'^ . 

However the proof of this facts is involved and, as such, beyond the scope of this 
presentation. Therefore we shall restrict ourselves to a succinct description of the 
principal steps towards the proofs. 

The basic idea is to show that if a metric M space admits a Wald curvature at any 
point, than M is locally homeomorphic to R^, thus any metric proprieties of 
can be translated to M, (in particular the first fundamental form). 
The first of these partial results is: 

Theorem 3.22. Let M be a convex metric space. Then M admits at most one 
Wald curvature Kwip), € M. 

Proof By Corollary 3.12. it suffices to prove that any disk neighborhood B{p ; p) e 
Af{p) contains a non degenerate sd-quad. 

Without loss of generality one can assume that B{p; p) contains three points pi , p2 , 
s.t. d{p,pi) < p/2, i — 1,2,3.^^ Then, by the convexity of M it follows that 
3g e M s.t. p ^ P2,V3 and p2q + p^q = P2P3- But P2P3 < PP2 + PP^ < P =^ 
(pq < p/2) V {pp2 < p/2). In the first inequality holds, then pq < pp2 +P2q < P, i-e 
q G B{p ; p); and if the second one holds, then pd < pps + p^q < p, i.e. q G B{p; p). 
But p ^ q, therefore p,P2,P3,q are not linear. 

□ 

Our next step will be to analyze those neighborhoods that display "a normal 
behavior" , both metrically and curvature- wise: that is precisely those disk neigh- 
borhoods in which the Wald curvature is defined and ranges over a small, bounded 
set of values prescribed by the very radius of the disk: 

Definition 3.23. A disk neighborhood B{p;p);p > is called regular iSV non- 
degenerate quadruple Q C B{p;p) , kw{Q) exists and |ki4'('9)| < 7r^/16p^. 



16 See ini 
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Remark 3.24. If Kwip) exists, then for any sufficiently small p, B{p]p) will be 
regular. 

It turns out that regular neighborhoods, in compact, convex spaces have the 
following "nice" (i.e. Real Plane like) proprieties: 

Proposition 3.25. Let M be a compact, convex metric space and let B{p;p) C 
M be a regular neighborhood. Then if a non- degenerate quadruple Q C B{p;p) 
contains two linear triples of points, then Q is linear. 

Proposition 3.26. Let M be a compact, convex metric space. Then: 
y p G M and ^ B{p ; p) regular, 3 g, r e B{p ; p) s.t. p, q, r are not linear. 

Proposition 3.27. Any regular neighborhood B(j); p) of a compact, convex metric 
space is strictly convex, i.e. q,r G B{p;p) ==> int{qr) C B{p;p). 

While the proof of this last Proposition is lengthy, that of the following important 
Corollary is not: 

Corollary 3.28. Let B(p; p) be a regular neighborhood. Then,y q,r € B{p; p), 3\ qr 
and int(qr) C B{p; p). 

Proof. By the convexity of B{p\ p) it follows the existence of at least one geodesic 
qr, y q,r e B{p;p). If s G int{q)r, then by the proposition above we have that 
s £ B{p ; p). It follows that B{p ; p) contains all the geodesies with end points q, r. 
Hence, by Proposition 3.25. , the geodesic segment qr is unique. □ 

We can now begin to prove that a compact, convex metric space locally mimics 
M^. We start by showing that the sinus function is defined on M, thus allowing for 
angle measure (hence for the definition of Polar Coordinates on regular neighbour- 
hoods^^). 

First, let M be as before, and let p € M s.t. kwP exists. Let q,r £ B[p ; p),q^ p ^ 
r, where B(j) ; p) is a regular neighborhood of p. Then, Vx S [0, min{pq,pr}), define 
q{x) S pq,r{x) £ pr by: d(p,q{x)) = x = d{p,r{x)), and let d{x) — d{,q{x),r{x)) 
(see Figure 10 bellow). 

Proposition 3.29. The following limit exists: 



x^O X 

We omit the proof since it is rather involved (but canonical for any axiomatic 
approach to Euclidian Geometry - see, for instance, |B], |RR| . ) 
Now we can define the measure of angles at p : 

Definition 3.30. The measure of the angle Z{q,p,r)) is given by: 



Remark 3.31. The Definition above enables us to define Polar Coordinates on reg- 
ular neighborhoods^* in the following manner: 

Let Pi,P2 G B[p;p) s.t. p,Pi,P2 a-rc not coUinear. (Such points exist by Propo- 
sition 3.26.). To every point q E B{p;p) we associate the following pair of real 

^^In the same way geodesic polar coordinates are used on classical surfaces. 

... and once Coordinates (be they Polar or Cartesian) are introduced, the (local) homomor- 
phism with is immediate. 
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Figure 10. 



numbers (defining the Polar Coordinates of q relative to the frame determined by 
P,Pi,P2)- (r{q),e{q)), where 



A 



r{q) = d{p, q) 



a( m{Z{q,p,pi)) i{\TniZ{p2,p,pi)) -m{Z{q,p,pi))\=miZiq,p,pi)) ; 



and 

^ 2tt - m{Z{q,p,pi)) if \m{Z{p2,p,pi)) - m{Z{q,p,pi))\ ^ Tn{Z{q,p,pi)) 

We can now safely state the foretold homomorphism result: 

Proposition 3.32. Any convex, compact metric space is locally homeomorphic to 
the real plane. 



4. Computing Embedding Curvature 

In this section we develop formulas for the computation of Embedding Curvature 
of Quadruples. First we follow the classical approach of Wald-Blumenthal that 
employs the so-called Cayley-Menger determinants (see bellow). Unfortunately, 
the formulas obtained, albeit precise are transcendental. Therefore we present, in 
the next subsection, the approximate formulas developed by C.V. Robinson. 



4.1. Embedding Curvature The Determinant Approach. Given a general 
metric quadruple Q ~ Q(pi,P2,P3,P4), of distances 

dij = dist{pi,pj), i = 1, ...,4, we denote by D{Q) — -D(pi,P2,P3,P4) the following 
determinant: 



(4.1) 



D{pi,P2,P3,Pi) 






1 


1 


1 


1 


1 





"12 


d'^ 
"l3 


dl 


1 


dh 





"23 


"24 


1 


dh 


"23 





dl. 


1 


dl. 


"24 


"34 
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Then the embedding curvature k{Q) of Q is given depending upon the embedding 
space (i.e. upon the sign of the curvature) - by the following formulae: 



(4.2) k{Q) 



if£)(Q)=0; 
K, K <0 if det{cosh k - dij) = ; 
K, K > if det{cos ^/n ■ dij) and y/n ■ dij < tt 

and all the principal minors of order 3 are > 0. 



The determinant D{Q) = -D(pi,P2,P3,P4) is called the Cayley-Menger determi- 
nant (of the points pi, ...^4)^^ and, in order to prove (4.2) we need first to investigate 
some of its properties. 
We start with the following 



Then: 
1 



Lemma 4.1. Let: pi, ...,p4 be points in ' 

(4.3) Gram{prp^,pip^,p^) = ^D{pi,p2,P3,P4) ; 

where 

Gram{piP2,PiP3,PiP4) = det{pipi ■ pipj)ij^2,3,4 



(Here "■" denotes the standard scalar (dot) product in R^.) 

Proof. Use expansion and manipulation of determinants. □ 

Since is a known fact that: 

Gram ipip2,piPa,pipi) = {Vol{pi,p2,p?„Pi)f ; 

where V^ol(pi,p2,P3,P4) denotes the (un-oriented) volume of the parallelepiped de- 
termined by the vertices pi, ...,P4 (and with edges P\P2,Pipi,p-Lpi)', formula (2.3) 
shows that: 

(4.4) D{pi,p2,P3,P4) = 8(Fo/(pi,P2,P3,P4))^ • 

Therefore the following assertion is immediate: 

Proposition 4.2. The points pi,...,p4 are the vertices of a simplex in M.^ iff 

D{pi,P2,P3,P4) ¥= 0: 

However, we can prove the much strong result bellow: 



Theorem 4.3. Let dij > , 1 < 4 , i ^ j. 

Then there exists a simplex T — T{pi, ...,^4) C 
DiPi,Pj) <0,{y){i,j} C {1,...,4} and 
D{pi,pj,pk)>0,{y) {i,j,A;} C {!,..., 4}; 
where, for instance, 



D{pi,P2) = 



s.t. dist{xi,Xj) = dij ,ij^j;iff 



1 


d% 



d^ 



12 




and 



Dipi,P2,P3) = 



1 dl2 



dl2 



^13 ^^23 



d?3 
"23 





This definition readily generalizes to any dimension, as do the results bellow. 
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etc... 

In fact, the necessary and sufficient condition above can be relaxed, indeed one 
can also show that the following holds^°: 

Proposition 4.4?^ Let rf.y > , 1 < 4 , i 7^ j. 

Then there exists a simplex T — T{pi, ...,^4) C R-^ s.t. dist{xi,Xj) = dij ,i ^ j; iff 
D{pi,p2,P3,P4) 7^ and sign D{pi,p2,P3,P4) = +1 • 

Proof. (Sketch) Sufficient to show (by using standard operations on determinants) 
that: 

D{pi,P2,P3)D{pi, ...,Pi, ...,P4) = Mil + D{pi, ...,p^, ...,P4)D{pi,p2,P3,P4) 

where Ma is the cofactor (in D) of and were we used the notation: 

{pi, ...,pi, ...,P4} = {P1,P2,P3,P4} \ {Pi}. □ 

Proving the formula for the spherical and hyperbolical cases would prove to be 
to technical for this limited exposition; suffice to say that they essentially reproduce 
the proof given in the Euclidian case, and tacking into account the fact that per- 
forming computations in the spherical (resp. hyperbolic) metric one has to replace 
the distances dij by cos dij (resp. coshd^)^^ 



4.2. Embedding Curvature Approximate Formulas. The formulas we just 
developed in are not only transcendental, but also the computed curvature may 
fail to be unique (see preceding section). However, uniqueness is guaranteed for 
sd-quads. Moreover, the relatively simple geometric setting of sd-quads allows for 
the development of simple (i.e. rational) formulas for the approximation of the 
Embedding Curvature. 

Proposition 4.5. Given the metric quadruple Q — Q{pi,p2,P3,P4), of distances 
dij = dist{pi,pj), i — 1,...,4, the embedding curvature k,{Q) is well approximated 
by: 

f4 5) K(0) = 6(cosZo2 + cosZo20 

^' d24(di2sin2(Zo2)+d23sin2(Zo2')) 
where: Zo2 — ^{piP2Pa) , Zo2' — Z(p3P2P4) represent the angles of the Euclidian 
triangles of sides di2, di4, d2i and d23, c?24, c?34 , respectively. 
The error R can be estimated by using the following inequality: 

(4.6) \R\ = \R{Q)\ = \k{Q) - K{Q)\ < AK^{Q)diam^{Q)/X{Q) 

where we put: A(Q) = 1^24(^12 sin Zo2+(i23 sin Zo2')/S''^, and where S — Max{p,p'}; 

2p = di2 + dl4 + ^24 , 2p' = ^32 + ^34 + ^24 • 

Proof The basic idea of the proof is to recreate, in a general metric setting, the 
Gauss Map - in this case one measures the curvature by the amount of " bending" 
one has to apply to a general planar quadruple so that it may be "straightened" 
(i.e. isometrically embedded as a sd-quad) in some S^. 

Consider two plane^^ triangles ApiP2P4 and AP2P3P4, and denote by ^Pi.kP2.kP4,k 



^'^For the direct proof of Theorem 4.3., see [Be) or, alternatively (b1 

^^We formulate this result - for convenience and practicality - for the case n = 3, only. 
However it is readily generalized to any dimension. 
22 See ini for the full details 
^'^i.e. embedded in B? = So 
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and Ap2,kP3,kP4:,k their respective isometric embeddings into Sk- Then pi^kPj,k 
wiU denote the geodesic (of Sk) through pi,k and pj^k- Also, let Zk2 and Zfe2' 
denote, respectively, the following angles of l^pi.kP2,kPA.k and Ap2.kP3.,kPi,k ■ ^fc2 ~ 
^Pi.kP2,kP4,k and Zfc2' = Zp2MP3,kP4,k- (See Fig. 11) 




Figure 11. 



But Zk2 and Zfc2' are strictly increasing as functions of k. Therefore the equation 

(4.7) Zfc2 + Zfc2' = vr 

has at most one solution k* , i.e. k* represents the unique value for which the points 
Pi,P2,P4 are on a geodesic in Sk (for instance on P1P4). 

But that means that k* is precisely the Embedding Curvature, i.e. k* = k{Q) , 
where Q = Q{pi,P2,P3,Pi)- 
Equation (4.7) is equivalent to 

2 '^k* 2 2 ^fe* 2 -, 

cos + cos — 1 

2 2 

The basic idea being the comparison between metric triangles with equal sides, 
embedded in Sq and Sk , respectively, it is natural to consider instead of the previous 
equation, the following: 

(4.8) e{k, 2) • cos^ ^ + e{k, 2') • cos^ ^ = 1 
where we denote: 

, , COS^^ , ,^ C0S2^ 

0{K2):= e{k.2'):^ 

cos^ cos"^ 

Since we want to approximate k{Q) by K{Q) we shall resort - naturally - to 
expansion into MacLaurin series. We are able to do this because of the existence 
of the following classical formulas: 

2 Zfc2 _ sin(p\/fc) ■ sin(d\/fc) 

COS ~ — ^ . — 5 /c U J 



2 sin(di2'\/^) • sin(d24'\/^) 
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2 Zfe2 smh{p\Ot) ■ sinh((i-\/fc) 

cos — — = j= ; K < ; 

2 sinh((ii2V A;) • sinh(rf24vfc) 



and, of course 

2 ^o2 pd 



cos 



2 di2d24 

were: d = p — di4 = {di2 + d24 — du)/'^ -^^ 

By using the development into series of f\{x) = ^'^^ and f2{x) = one 
(easily) gets the desired expansion for 6{k,2): 

(4.9) e{k, 2) = 1 + ^ fcrfi2rf24 ( cos(Zo2) - 1) + r ; 

where: |r| < |A:^p'' , for \kp^\ < 1/16. By applying (4.9.) to (4.8), we receive: 

(4.10) [1 + ^rdi2rf24(cos(Zo2) - 1) +r]cos2 ^ + 

[1 + ifc*d23d24(cos(Zo2') - 1) + r'] cos2 ^ = 1 ; 
for: \r\ + \r'\ < l{k*)^{Max{p,p'})^ = ^ik*fS^ . 

By solving linear equation (in variable k*) (4.10) and using some elementary trigono- 
metric transformation one has: 

^* ^ 6(cosZo2 + cosZo20 ^ ^ 

rf24(rfi2 sin2(Zo2) + d23 sin2(Zo2')) 

where: 

,p, ^ 12(|r| + |/|) 9(fc*)2max{p,p'} 



rf24(di2 sin^(Zo2) + das sin2(Zo2')) ^24(^12 sin^(Zo2) + ^23 sin2(Zo2')) 
But k* = k{Q) so we get the desired formula (4.5) . 

To prove the correctness of the bound (4.6) one has only to observe that: 
S = Max{p,p'} < 2diam{Q), (diam{Q) = max {rfj,}), 

l<i<j<4 

and perform the necessary arithmetic manipulations. 

□ 

Remark 4.6. (a) The function A = X{Q) is continuous and 0-homogenous as a 
function of the dij-s. Moreover: \{Q) > and X{Q) = ^ sinZo2 = sinZo2' = 0, 
i.e. iff Q is linear. [Therefore for sd-quads X(Q) > and, moreover, X{Q) ^ =J> 
Q linearity.] 

(b) Since X{Q) 7^ it follows that: K{Q) e M for any quadrangle Q. 
In addition: sign(k(Q)) — sign(K{Q)). 

(c) If Q is any sd-quad, then {Q)diam'^ {Q) / X{Q) < 00. Moreover, 
if X{Q) :^ 0,2^ then k'^ {Q)diam'^ {Q) / X{Q) ^ i.e. \R\ is small if Q is not close to 
linearity. 

In this case |i?((3)| ^ diam?{Q) (for any given Q). 



■^^and the analogous formulas for cos^ ~2~~' 
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i.e. for not very small values of 
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Since the Gaussian curvature kaip) at a point p is given by: 

kaip) = lim K{Qn) ; 

where Qn Q = Opipp^p^ ; diam{Qn) O,from Remark 4-6.(c) we immediately 
infer that the following holds^^: 

Theorem 4.7. Let S be a differentiable surface. Then, for any point p G S: 

kaip) = lim K{Q„) ■ 

rx— »0 

for any sequence {Qn} of sd-quads that satisfy the following condition: 
Qn^ Q = Opipp3P4 ; diam{Qn) . 

Remark 4.8. In the following special cases even "nicer" formulas are obtained: 

(1) If rfi2 = ^32, then 

^ ' ^ ' di3 ■ d24 sin^ Zo2 + sin^ Zo2' 

(here we have of course: dis = 2dx2 = 2^32); or, expressed as a function of 
distances alone: 

(A m K(0) - 12 2df2 + 2d^4 - - 

^ ' &dl,dl, - [dj, + dl, dl,Y {dl, + dl, - dl,f 

(2) If d\2 = d32 = d24 and if the following condition also holds: 

(3) Zo2' = 7r/2; i.e. if ^34 = + <^24 or, considering (2), also: ^34 = 2(if2 
then 

f4 13) K(0) - - ^^12 - ^14 

^ ' di2(l + sinVo2) Mi2+Ml,dl^-df,- 

5. Appeendix 1 - The Menger and Haantjes Curvatures 

Better known than the Wald Curvature, the Menger Curvature is a metric def- 
inition of curvature of curves, as is the Haantjes Curvature. As such they can be 
employed as sectional curvatures to approximate curvature of triangulated surfaces. 
We begin by introducing the the Menger Curvature: this is a metric expression for 
the circum-radius of a triangle^^, based upon elementary high-school formulas: 

Definition 5.1. Let (M, d) be a metric space, and let p,q,r G M he three distinct 
points. Then: 

„ , . \/ (pq + qr + rp){pq + qr - rp){pq - qr + rp){-pq + qr + rp) 

^M(P,Q,r) = ; 

pq - qr ■ rp 

is called the Menger Curvature of the points p, q, r. 

We can now define the Monger Curvature at a given point by passing to the 
limit: 

Definition 5.2. Let (M,d) be a metric space and let p € M be an accumulation 

point. Then M has at p Menger Curvature km{p) iff 

Ve > 0, 36 >0 s.t. d{p,p^) <S; i = 1,2,3 =^ \K{Q) - km{p)\ < e. 



^''and gives theoretical justification to the algorithm 

^'^ thus giving in the limit a metric definition of the Osculatory Circle 
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Remark 5.3. The apparent equivalent notion of Alt Curvature, in which one uses 
only two points converging to the third, is in fact more general, where we define 
the Arp curvature by: 

Definition 5.4. Let (M,d) be a metric space and let P £ M be an accumulation 
point. Then M has at p Alt Curvature ka{p) iff the following limit exists 

KAip) = lim K{p,q,r) . 

q.r^p 

However, both km{p) and ka{p) suffer from the same imperfection: since they 
are both modelled closely after the Euclidian Plane, they convey this Euclidian 
type of curvature upon the space they are defined on. However, the next definition 
doesn't mimic closely so it better fitted for generalizations: 

Definition 5.5. Let (M,d) be a metric space and let c : / = [0, 1] ^ M be a 

homeomorphism, and let p,q,r g c(/), q,r ^ p. Denote by qr the arc of c(/) 
between q and r, and by qr segment from q to r. (See Figure 12 bellow.) 




Figure 12. 

Then c has Haantjes Curvature kh{p) at the point p iff: 
.lip) = 24 lim ; 

where "l{qr)" denotes the length^* of qr. 

Remark 5.6. kh exists only for rectifiable curves, but if km exists at any point p 
of c, then c is rectifiable. 

Remark 5.7. Evidently we have the following relationship between curvatures: 

3 KM =^ 3 KA ■ 

while 

3 KA 3 KM ■ 

However, we can prove the following theorem: 



given by the intrinsic metric induced by d 
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Theorem 5.8. Let c : I ~* AI be a rectifiable curve, and let p G M . 
If KA ( or Km ) exists, then kh (p) exists and 

Remark 5.9. This last result and the Remark preceding it allow as to employ any 
of the curvatures above in estimating curvatures of smooth curves on triangulated 
surfaces. 

6. Appendix 2 - The Rinow Curvature 

The curvatures introduced before may seem a bit archaic in comparison to the 
more fashionable approach of comparison triangles, with their ar reaching applica- 
tions. We present here one of these comparison criteria and show its equivalence 
with the Wald curvature. We start with the following definition: 

Definition 6.1. Let {M,d) be a metric space, together with the intrinsic metric 
induced by d. Let R — int{R) C M be a region of M . We say that i? is a region of 
curvature < k (k G R) iff 

(1) Vp, (7 e i? , 3 a geodesic segment pq C R; 

(2) yT{p, q,r) C R is isometrically embeddable in 5^; 

(3) If T{p, q,r) C R and x &pq,y G pr, and if the points Pi^,qn,r^,Xf,,yi^ e S,^ 
satisfy the following conditions: 

(a) T{p,q,r)^T{p,,q^,r,); 

(b) T(p,q,x) = T{p,„qf„x^); 

(c) T(p,r,y) ^T{p^,r^,yi^); 
then xy < x^y^. 

By replacing the condition: xy < XK,yK^ with: "a;?/ > x^J/k", we obtain the defini- 
tion of a region of curvature > k. (See Fig. 13.) 




Figure 13. 
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We now pass to the localization of the Definition above: 

Definition 6.2. Let (M, d) be a metric space, together with the intrinsic metric 
induced by d, and let p G M be an accumulation point. Then M has at p Rinow 
Curvature hr{p) iff 

(i) eAf{p), linear; 

(ii) Ve > 0, 3(5.0, s.t. B{p; 6) is (a) a region of Rinow curvature < kr{p) + s and 
(b) a region of Rinow curvature > k,b.{p) — e. 

While its greater generality endows the Rinow curvature with more flexibility 
in applications and makes it easier in generalization, it is even more difficult to 
compute than Wald Curvature. However this quandary was has an almost ideal 
solution, due to Kirk (see |K1)) solution which we briefly expose here: 

Definition 6.3. Let Af be a compact, convex metric space, and let p S M . 
ff K,w{p) exists, then Kji{p) exists, and kr(j>) = Kwip)- 

Unfortunately, since kb,{p) makes no presumption of dimensionality, the exis- 
tence of kr{p) does not imply the existence of Kwip)- 

Counterexample 6.4. Let M = M^. Then kr{p) = but kw{p) does not exist at 
any point, since every neighborhood contains linear quadruples. 

The solution (due to Kirk) of this problem is to consider the Modified Wald 
curvature nwK, defined as follows: 

Definition 6.5. Let (M, d) be a metric space, together with the intrinsic metric 
induced by d, and let p Cz M . Then M has at p Modified Wald curvature Curvature 
i^wk{p) iff 

ii)^N e Af{p), N hnear; 

(ii) Ve > 0, 3(5.0, s.t. if Q C B{p;S) is a non-degenerate sd-quad, then K^r{Q) 
exists and \kwk{p) — kw{Q)\ < £■ 

Remark Q. 6. 3kw{p) 3 kwk{p) ^'^t ^kkw{p) kw{p)- 

This modified curvature indeed represents the wished for solution, as proved by 
the following to Theorems: 

Theorem 6.7. Let {M,d) be a metric space. Then: 

3kr{p) =^ 3kwk{p) and kr{p) = kwk{p)- 

Theorem 6.8. Let (M, d) be a metric space together with the associated intrinsic 
metric, and let p M . Then, if 

(i) kwk{p); 
and if 

(ii) 3B{p-p) e N{p), s.t. qr C B{p;p), \/q,r e B{p;p); 
then kr{p) exists and kr{p) = kwk{p)- 

7. Appendix 3 - The Radius Formula 

The Cayley-Menger determinant allows one to express not only the volume and 
area^^ of simplices in M" but (as expected) it may be used to compute the radius 

^^The 2-(iimonsional analogue of Formula (4.4) for the area of the triangle T(pi,p2iP3) being: 
[Area (pi,P2,P3))^ = -D{pi,p2,P3) ■ 
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of the circumscribed sphere around an Euchdian simplex. To be more precise, we 
have the following result^*^: 

Theorem 7.1. (1) The radius R = R{pi,p2,P3,P4) of the sphere circumscribed 
around the tetrahedron T{pi,p2,p3,p4) € is given by: 

1 A{pi,P2,Pz,pa) 



R^ 



2D{pi,p2,P3,P4) 



where: 



A(pi,P2,P3,P4) 



(2) The points Pi,P2,P3,P4,P5 S are coplanar or co-spherical iff 






dl2 


rf?3 


dl4 


dl2 
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"24 


dl3 


"23 





dl4 


dl4 


"24 
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Proof 



(1) If e 

obtain: 



A(pi,P2,P3,P4,P5) = . 

is s.t. rfoi = R, i = 1,...,5, then by direct computation we 



r(po, •••,P5) = ~2R^T{pi, ...,P5) - A(pi,P2,P3,P4,P5) ; 

from which the desired formula follows immediately if we chose po as the 
center of the sphere circumscribed around the points Pi,P2,P3,P4,P5- 
(2) (^^) Let {x^,x^,x^} be any orthonormal coordinate frame for M?, and 
lot pj: i = 1, .... 5: j = 1, 2, 3: represent the coordinates of the points Pi,P2,P3,P4^P5 
relative to this coordinate system. Then Pi,P2,P3,P4,P5 belong to the same sphere 
or plane iff 3 (a, b, cj) ^ (0, 0, 0), j = 1, 2, 3; s.t. 



aWPiW^ +b + ^Cjp'- ; i = 1,...,5. 
Then Ai = A2 = 0, where: 

Al(pi,P2,P3,P4,P5) 



and 



^2{Pl,P2,P3,P4,P5) 





1 p\ 


P\ 


Pi 


\\P2\\' 


1 pI 


pI 


pI 




1 pI 


pI 


pI 


\\P4\\' 


1 pI 


p\ 


pI 


\\P^\? 


1 Pi 


p\ 


pI 


IbilP 


-2pl 


-Ipl 




\\P2\\' 


-M 


-2pl 




\\P2\? 


-'2pI 


-'^pI 




\\P2\\' 


-M 


-2pl 




\\P2\? 




-2pl 





2pl 



2pi 



Therefore Ai • A^ = 0. But Ai • A^ = A(pi,p2,P3,P4,P5), so this implication is 
proven. 



^'^that can be readily generalized to higher dimensions 
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(<J=) A(pi,p2,P3,P4,P5) = =4> Ai = and there exist numbers {a,b,Cj) ^ 
(0,0,0), i = 1,2,3; s.t. 

3 

a||Piir+^ + I]cjP^ =0; i = l,...,5; 

i=i 

i.e. Pi,P2jP3,P4,P5 belong to the plane or the sphere given by the equation 

3 

a\\X\f + b + cY,CjX = 0; X = {xuX2,X3) . 

□ 
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